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摘 要 :大 型 飞机 结构 弹性 会 显著 的 影响 到 机 体 的 载荷 大 小 及 分 布 ,准确 计算 结构 弹性 引起 的 载荷 
变化 对 提高 飞机 设计 水 平 具有 重要 意义 。 首 先 给 出 了 一 种 机 翼 气 动 载荷 静 气 弹 理 论 修 正方 法 ,该 
方法 应 用 工程 梁 理 论 计 算 机 杏 变 形 , 应 用 涡 格 法 计算 变形 引起 的 气动 力 增 量 ,然后 进行 迭代 直至 收 
敛 , 用 升力 面 理论 对 NACA -TN3030 报告 中 的 升力 线 理论 进行 了 改进 提高 。 通 过 风 洞 试验 证 明了 
该 理论 方法 的 准确 性 。 在 飞行 试验 中 测量 得 到 了 机 杏 变 形 和 飞行 参数 两 方面 的 数据 ,分 别 以 这 两 
方面 的 数据 为 输入 各 自 进行 弹性 修正 计算 ,两 种 途径 得 到 的 结果 相互 验证 ,修正 结果 吻合 较 好 ,经 
过 载荷 修正 后 翼 根 谊 答 降 低 了 大 约 3% ,该 方法 易于 在 工程 中 推广 使 用 ,对 于 减轻 机 杏 结 构 重 量 有 


显著 意义 。 
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Abstract :The structural elasticity of large aircraft will significantly affect the load size and distribution of 
the airframe. lt is important to calculate the load variation caused by structural elasticity accurately for 
improving the design level of aircraft. In this paper，a modified method of aeroelastic theory for 
aerodynamic load of wing is presented, which uses the theory of engineering beam to calculate the 
aerodynamic increment caused by the deformation , then iterates to converge , and improves the lift line 
theory in NACA - TN3030 report by using the theory of lift surface. The accuracy of the method is proved 
by wind tunnel test. In flight test,the data of wing deformation and flight parameters are measured ,and the 
elastic correction calculation is carried out respectively with the data of these two aspects as input. The 
results of the two methods are verified by each other,and the correction results are in good agreement with 
each other. After load correction ,the bending moment of wing root is reduced by about 3% . This method is 
easy to be popularized in engineering. lt is significant to lighten the weight of wing structure. 
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一 架 飞 机 从 起 飞 到 着 陆 全 过 程 无 时 无 刻 不 在 承 
受 着 载 集 ,包括 空气 动力 载荷 发 动机 推力 和 上 自 刁 的 
惯性 载荷 。 载 荷 是 静 强 度 、 动 强度 和 效劳 强度 分 析 
的 重要 输入 ,准确 确定 结构 载荷 是 飞机 结构 强度 设 
计 的 基础 。 如 果 设 计 载 荷 小 于 实际 使 用 情况 ,必然 
因 结 构 破 坏 导 致 飞行 事故 ;反之 否 设 计 载 傈 过 大 则 
会 导致 飞机 结构 重量 增加 而 影响 飞机 的 性 能 。 在 我 
国 大 飞机 工程 发 展 之 初 ,为 保证 我 国 研制 的 大 型 飞 
机 具有 国际 竞争 力 ,国内 专家 就 将 大 型 飞机 载 谷 确 
定 技术 列 为 必须 突破 的 关键 技术 之 一 。 

随 着 飞机 设计 技术 的 不 断 发 展 ,飞机 性 能 不 断 
提高 。 对 于 大 型 运输 类 飞机 ,为 了 提高 经 济 性 ,出 于 
增 升 减 阻 的 目的 ,通常 采用 更 大 展 弦 比 的 机 嗓 较 薄 
的 辟 型 ,为 了 增加 座位 而 加 长 机 屿 ,这 些 使 得 飞机 变 
得 更 长 .更 细 、 更 柔软 ”1。 因 此 ,飞机 的 弹性 变形 对 
载荷 的 影响 也 越 来 越 重要 ,准确 地 计算 弹性 飞机 在 
实际 飞行 中 的 载荷 ,对 于 飞机 的 结构 和 强度 设计 非 
党 重要 。 大 型 飞机 的 机 如 在 飞行 过 程 中 会 产生 显著 
的 弹性 变形 ,此 变形 会 对 机 嗓 的 气动 载 答 分 布 产 生 
明显 的 影响 。 在 国 军 标 ” 和 适 航 规范 ”中 都 有 相 
应 的 要 求 ,要求 在 载荷 计算 中 必须 考虑 结构 变形 对 
载荷 分 布 的 影响 ,也 就 是 要 进行 弹性 修正 。 另 一 方 
面 ,现代 大 型 运输 机 通常 采用 大 展 弦 比 的 后 掠 机 豆 ， 
该 类 型 的 机 如 载 傈 经 过 弹性 修正 后 会 有 较为 明显 的 
降低 ,这 对 于 减轻 机 缀 结构 重量 很 有 帮助 。 

按照 规范 要 求 , 飞 机 结构 载 答 的 种 类 众多 ,经 过 
多 个 飞机 型 号 研制 的 实践 证 明 , 稳 定 俯仰 机 动 是 机 
让 载荷 的 设计 情况 ,因为 在 该 机 动 中 全 机 不 仅 要 达 
到 最 大 的 法 向 过 载 ,同时 平 尾 上 又 有 很 大 的 向 下 的 
配 平 载 知 。 因 此 对 稳定 俯仰 情况 的 机 避 载 谷 进 行 修 
正 后 就 可 以 降低 机 辟 的 设计 载荷 ,而 其 他 类 型 的 机 
动 中 ,机 距 载 荷 显 著 小 于 稳定 俯仰 机 动 ,进行 修正 的 
意义 不 大 。 稳 定 俯仰 机 动 是 对 称 机 动 的 一 种 ,飞机 
在 1g 平 飞 的 初始 状态 下 以 很 缓慢 的 速率 偏转 升降 
舵 ,俯仰 角 加 速度 为 零 。 在 此 机 动 中 ,飞机 总 是 处 于 
平衡 状态 ,变形 对 于 结构 固有 频率 非常 缓慢 ,可 以 忽 
略 变形 引起 的 惯性 力 ,因而 该 问题 可 以 用 静 气 弹 的 
理论 进行 分 析 。 

从 很 早 之 前 人 们 就 开始 研究 机 可 载 荷 的 静 气 弹 
修正 方法 。 近 些 年 随 着 计算 机 技术 的 发 展 ,对 于 静 
气 弹 问题 越 来 越 多 的 学 者 采用 采用 计算 结构 力学 
(computational structural dynamics ,CSD ) 与 计算 流体 
力学 (computational fluid dynamics, CFD ) 耦合 计算 的 
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方法 研究 机 杜 的 静 气 动弹 性 响应 问题 。 计 算 中 通常 
采用 非 结 构 动 网 格 技术 来 适应 迭代 过 程 中 的 机 骂 变 
形 '“ ,通过 RBF( 径 向 基 范 数 ) 方 法 实现 流 固 耦合 
面 的 数据 交换 '" 1 。 在 描述 结构 变形 时 ,也 有 采用 模 
态 到 加 的 方法 “; 。CFDZCSD 耦合 的 方法 需要 较 大 
的 计算 耗费 ,不 利于 工程 应 用 ,工程 上 目前 更 多 的 仍 
是 采用 相对 简化 但 精度 满足 工程 要 求 的 做 法 。 早 年 
的 NACA -TN3030 给 出 来 一 种 适合 于 工程 应 用 的 
修正 方法 ” ,其 核心 思想 是 :在 刚性 载荷 的 基础 上 ， 
利用 工程 染 理 论 计 算 机 辟 的 变形 ,用 升力 线 理论 计 
算 变 形 引起 的 气动 力 增 量 ,然后 在 这 二 者 之 间 进 行 
迭代 直至 收敛 。 这 种 方法 具有 物理 概念 清晰 、 计 算 
快捷 等 优点 。 国 内 在 工程 实践 中 ,采用 这 种 方法 研 
究 了 静 气 弹 效 应 对 机 翼 载 荷 大 小 及 分 布 的 影 
响 W "1, 但 是 这 种 方法 的 气动 力 计算 采用 的 升力 线 
理论 对 于 较 大 后 掠 角 的 机 可 会 有 比较 大 的 误差 。 同 
时 对 于 静 气 弹 的 研究 更 多 关注 的 是 其 于 对 气动 和 飞 
行 性 能 的 影响 ,对 载荷 的 关注 较 少 ,关于 静 气 弹 载 荷 
分 布 修正 的 风 洞 试验 很 少 '" ,通过 飞行 试验 进行 研 
究 的 更 是 还 没有 公开 的 发 表 的 文献 。 

本 研究 首先 给 出 了 一 种 理论 分 析 方 法 ,使 用 升 
力 面 理论 计算 气动 力 增 量 , 从 而 对 NACA TN3030 中 
的 方法 进行 了 改进 ,其 次 通过 风 洞 试验 证 明了 理论 
分 析 方 法 的 准确 性 ,最 后 使 用 在 试飞 试验 中 得 到 的 
两 方面 实测 数据 。 应 用 该 理论 分 析 方 法 通过 两 种 途 
径 对 机 可 气 动 载荷 进行 了 静 气 弹 人 和 修正, 两 种 途径 的 
修正 结果 可 以 相互 比较 验证 。 


1 理论 分 析 


静 气 弹 计算 通常 需要 进行 结构 变形 和 变形 引起 
的 气动 力 增 量 计算 ,并 且 要 在 两 者 之 间 进行 迭代 直 
至 收敛 ,下 面 分 别 叙 述 。 


1.1 机 翼 变 形 计算 原理 


对 于 大 展 弦 比 的 机 既 , 可 以 用 一 根 变 刚度 的 工 
程 梁 来 描述 “ 。 梁 的 抗 弯 刚度 和 抗 扭 刚度 分 别 为 
E1(z) 和 GJ(z) ,z 坐标 是 沿 机 嗓 刚 轴 由 忆 根 指向 疲 
兴 。 机 缀 的 变形 用 梁 的 变形 来 代替 ,机 姻 的 变形 有 
弯曲 变形 和 扭转 变形 。 根 据 莫 尔 定理 ,在 弹性 刚 轴 
上 任意 一 点 的 弯曲 转角 和 绕 刚 轴 的 扭转 角 可 以 用 单 
位 载 集 法 来 求 , 即 在 4 点 加 一 单位 载 集 ( 弯 矩 或 扭 
和 矩 ) ,然后 根据 其 引起 刚 轴 上 各 切面 的 载荷 分 布 和 刚 


投稿 网 站 :http://ejam. xjtu. edu. cn 微 信 公众 号 :应 用 力学 学 报 


第 6 期 


轴 实 际 所 受 的 载荷 即 可 求 得 弯曲 变形 角 和 扭转 变形 
角 ,其 表达 式 分 别 为 


_ FM (z)M(z) 
p4 = | Gy 
T(z) T(z) 
0 = | 本 | 


式 中 : gp 为 在 载荷 M(z) 作用 下 刚 轴 上 4 点 的 弯曲 
角度 ; 04 为 在 载荷 7(z) 作用 下 刚 轴 上 4 点 的 扭转 
角度 ; M”(z) 为 在 4 点 作用 一 单位 弯 矩 后 刚 轴 上 各 
截面 的 弯 矩 分 布 ; T(z) 为 在 4 点 作用 一 单位 扭矩 
后 刚 轴 上 各 截面 的 扭矩 分 布 ; £1(z) 为 机 辟 刚 轴 的 
弯曲 刚度 ; GJ(z) 为 机 经 刚 轴 的 扭转 刚度 。 

对 于 大 展 弦 比 的 后 掠 机 翼 由 于 机 翼 的 这 曲 和 捏 
转变 形 使 得 机 豆 的 顺 气流 剖面 的 迎 角 发 生变 化 。 对 
于 后 掠 机 可 , 符 刚 轴 的 后 掠 角 为 A, 机 波 弯 曲 角度 和 
扭转 角度 与 顺 气流 剖面 的 迎 角 改变 量 a, 之 间 的 关 
系 为 


Qa, = OcosA - psinA (3) 
式 (3) 表 明 , 机 可 谢 面 的 迎 角 改 变量 可 以 表示 
成 弯曲 变形 和 扭转 变形 的 函数 ,而 这 两 个 角度 又 与 
刚 轴 上 的 载荷 分 布 有 密切 的 关系 。 由 此 关系 可 以 建 
立 气动 弹性 变形 引起 的 迎 角 增 量 和 升力 分 布 之 间 的 
关系 。 
在 实际 的 计算 中 ,按照 波 肋 的 布置 ,将 机 灵 沿 展 
向 划分 为 若干 个 站 位 ,并 从 辟 根 到 翼 尖 顺序 编号 ,第 
i 个 站 位 的 弯曲 刚度 和 扭转 刚度 分 别 为 B1(i)、 
GJ(i) ,这 个 站 位 上 作用 的 弯 矩 和 扭矩 分 别 为 M(i) 
和 7T(i) , 则 式 (1)、 (2) 的 连续 积分 可 以 用 离散 的 代 
数 和 代替 ， 并 将 其 代入 式 (3) 中 得 
T(k)da(k M(k)d(k 
“AD 
式 中 : d(%) 为 沿 刚 轴 方向 两 个 肋 站 位 之 间 的 刚 轴 长 
度 。 在 实际 的 计算 过 程 中 , 净 载 和 荷 的 分 布 是 已 知 的 ， 
也 就 是 力 的 大 小 和 作用 点 是 已 知 的 ,结合 机 可 肋 站 
位 的 布置 和 刚 轴 的 位 置 ,就 可 以 确定 出 每 个 站 位 的 
弯 矩 和 扭矩 ,代入 式 (4) 即 可 得 到 每 个 肋 站 位 沿 来 
流 方 向 的 扭转 角 , 即 顺 气流 方向 迎 角 变 化 量 。 


1.2 结构 变形 引起 的 气动 力 计 算 


对 载荷 进行 静 气 弹 修正 时 气动 力 计算 的 是 变形 
引起 的 气动 力 增 量 。 本 研究 采用 涡 格 法 计算 , 涡 格 
法 是 升力 面 理论 中 一 种 比较 实用 的 计算 方法 ,目前 
在 气 弹 和 载荷 分 析 、 操 稳 、 飞 行 力学 和 飞 控 等 专业 中 
仍 被 广泛 使 用 "'“”l。 对 于 汉 身 组 合体 ,可 将 它 投影 


inA (4) 
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在 坐标 平面 上 按 一 个 当量 机 可 平 面 处 理 。 涡 格 法 在 
许多 文献 中 均 有 介绍 '” ,该 方法 比 NACA-TN3030 
中 使 用 的 升力 线 理论 在 计算 较 大 后 护 角 和 较 小 展 弦 
比 的 机 辟 载 荷 时 具有 更 高 的 精度 。 同 时 ,在 飞机 研 
制 的 详细 阶段 ,都 会 有 比较 完备 的 刚性 气动 力 数据 ， 
对 载荷 进行 弹性 修正 就 是 以 刚性 载荷 为 输入 ,计算 
出 结构 变形 引起 的 载 集 增 量 ,而 该 增 量 与 疲 型 厚度 、 
村 度 的 关系 并 不 大 ,因此 涡 格 法 应 用 在 本 方法 中 是 
可 行 的 。 

第 j 个 涡 格 上 的 马蹄 涡 在 第 i 个 涡 格 控制 点 处 
产生 的 y 向 无 量 纲 诱导 速度 值 是 习 ,气动 力 影响 系 


DV, 


数 C; 将 也 与 第 j 个 涡 格 上 的 渴 强 联系 起 来 ,由 毕 
奥 - 萨 瓦 公式 和 几何 关系 有 


Cs = (9 
7 


ov” y=1 


1 1 
re 
医 一 Xi ) (Xi — Xj) + (zz — 21) (2 — 21) 
V(x 一 xD) + (Z; — 21)” 
(xz 一 XU) (Xi — Xz;) + (zz 一 57) (2; | i 


V(x x )” + (2 -27)” 


1 Xi — XN 
?1 a 2 V(x 一 xD) + (Zi | 
1 | Xi 一 Ny 
221 Mo 一 和) 十 ( 关 2 | 加 
证 = (6) 
将 辟 面 线 化 边界 条 件 式 应 用 于 第 ;个 控制 点 上 
= [E02] -a (7) 


o0 


对 于 第 i 个 涡 格 ,根据 它 所 处 的 站 位 ,由 式 (4) 
得 到 对 应 的 a; , 计 人 所 有 训 格 对 它 的 影响 ,可 得 


>cm = Qi (8 ) 
求解 式 (8 ) 即 可 得 到 每 个 涡 格 上 的 涡 强 ,再 由 
茹 科 夫 斯 基 升 力 定 理 即 可 得 到 该 涡 格 上 的 气动 力 


增 量 


EO 


1.3 和 迭代 计算 


以 刚性 机 如 的 载 谷 作为 输入 ,计算 载荷 作用 下 
的 机 可 变形 ,然后 计算 此 变形 引起 的 气动 载荷 增 量 ， 
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将 此 增 量 锥 加 在 刚性 机 既 的 载荷 之 上 进行 下 一 次 达 
代 , 直 至 收敛 ,迭代 收敛 的 条 件 是 相 邻 两 次 的 迭代 各 
肋 剖 面 产 生 的 迎 角 增 量 不 再 变化 。 实 际 的 计算 过 程 
中 可 以 取 最 外 侧 辟 肋 位 置 的 顺 气 流 迎 角 不 再 变化 。 


1.4 理论 分 析 结 果 风 洞 试验 验证 


为 了 验证 上 述 理论 修正 方法 的 正确 性 ,加 工 了 
3 副 刚度 各 不 相同 的 机 可 进行 风 洞 试验 。 其 中 一 副 
是 刚性 机 避 , 另 外 两 副 是 弹性 机 桔 。 弹 性 模型 1 的 
刚度 由 实际 机 可 刚度 缩 比 得 到 。 弹 性 模型 2 的 刚度 
取 为 模型 1 的 70% ,通过 降低 模型 刚度 .增加 模型 的 
变形 量 来 进一步 验证 本 方法 。 在 风 洞 试验 中 ,对 各 
个 模型 进行 测 力 和 测 压 , 得 到 其 在 各 工 况 下 的 气动 
载 集 分布。 然后 以 刚性 模型 的 气动 载 全 和 弹性 模型 
的 刚度 数据 为 输入 ,按照 前 面 叙述 的 方法 对 刚性 载 
位 进 行 修正 ,得 到 弹性 模型 的 预计 的 载 傈 ,并 将 其 与 
试验 测量 得 到 的 实测 弹性 载荷 进行 对 比 , 即 可 验证 
方法 是 否 准确 。 图 1 和 图 2 分 别 给 出 了 两 副 模型 的 
对 比 结果 。 
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图 1 弹性 模型 1 的 弯 矩 理论 分 析 与 试验 实测 对 比 


Comparison between theoretical analysis and experimental 


Fig. 1 
measurement of bending moment of elastic model 1 
16 Ma=0.4, wa=10? 
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图 2 弹性 模型 2 的 弯 矩 理论 分 析 与 试验 实测 对 比 


Fig.2 Comparison between theoretical analysis and experimental 


measurement of bending moment of elastic model 2 
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同 马赫 数 和 迎 角 工 况 下 ,理论 分 析 得 到 的 弹性 模型 
的 弯 和 矩 数值 和 试验 实测 得 到 数值 吻合 程度 较 好 , 理 
论 计算 和 试验 实测 得 到 的 载荷 变化 量 的 趋势 和 数值 
吻合 较 好 。 证 明了 该 分 析 方 法 可 以 有 效 解决 工程 实 
际 中 的 载荷 静 气 弹 修正 问题 。 实 际 工程 应 用 中 有 两 
点 是 需要 注意 的 :中 气动 力 如 果 是 在 巡航 (也 就 是 
lg 平 飞 ) 状 态 下 给 出 的 , 则 需要 在 基础 载 集 中 减 去 
1g 平 飞 的 载荷 ,然后 再 进行 修正 ,否则 会 导致 修正 
量 偏 大 ;@) 修 正 后 机 枝 上 的 总 升力 会 有 所 降低 ,这 样 
就 不 满足 规定 的 过 载 系数 ,这 时 需要 将 整个 飞机 的 
迎 角 增 大 一 定 的 量 , 使 之 满足 过 载 系数 要 求 , 同 时 在 
机 可 载荷 上 需要 加 上 因 迎 角 增 大 引起 的 气动 载荷 


增 量 。 


2 试飞 验证 


2.1 试飞 试验 概况 


为 了 验证 本 研究 修正 方法 的 准确 性 ,规划 了 相 
应 的 飞行 试验 。 试 飞 飞 机 为 茶 大 型 飞机 ,机 中 外 形 
为 大 展 弦 比 后 掠 机 蛋 , 飞行 M, =0.6, H = 7000 m。 
在 飞行 试验 中 ,飞机 按照 规范 要 求 进行 各 种 机 动 , 在 
机 动 飞行 的 过 程 中 通过 相机 实时 测量 机 骂 上 辱 干 数 
量 足 以 描述 机 嗓 变 形 的 特征 点 的 坐标 ,通过 对 不 同 
时 刻 特 征 点 的 坐标 进行 比 对 ,就 可 以 得 到 不 同时 刻 
的 相对 变形 。 同 时 对 迎 角 、 侧 滑 角 、 控 制 面 的 偏转 角 
度 和 飞机 的 角速度 进行 测量 ,这 些 参 数 是 计算 气动 
载荷 的 关键 参数 。 图 3 给 出 了 右 机 史上 各 测量 点 的 
位 置 ,图 4 中 给 出 了 法 向 过 载 的 时 间 历 程 。 
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图 3 右 机 辟 位 移 测量 点 位 置 
Fig.3 Right wing displacement measurement point position 
对 于 大 型 飞机 ,在 飞机 设计 时 机 站 的 气动 载 从 
是 在 巡航 外 形 下 给 出 的 ,也 就 是 说 飞机 设计 时 气动 
载荷 的 计算 已 经 考虑 到 了 在 1g 过 载 时 的 机 缀 变形 ， 
因此 只 需要 对 在 1g 基础 上 的 增 量 载荷 进行 修正 ,机 
人 到 变形 量 也 相应 取 在 1g 基础 上 的 增 量 。 所 以 在 机 


由 图 1~2 可 以 看 出 ,弹性 模型 1 和 模型 2 在 不 


动 过 程 中 选取 了 2 个 典型 时 刻 点 ,其 中 A 时 刻 的 法 
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向 过 载 为 1.0g,B 时 刻 的 法 向 过 载 为 2.33g, 是 整个 
飞行 过 程 中 的 最 大 值 ,如 图 4 所 示 。2 个 工 况 之 间 
飞机 燃油 消耗 的 重量 相对 于 飞机 重量 来 说 是 一 个 很 
小 的 量 ,计算 中 忽略 不 计 。 图 5 给 出 了 机 票 上 各 位 
移 测量 点 在 B 时 刻 相对 A 时 刻 的 变形 量 。 表 1 给 
出 了 机 正典 型 部 位 的 相对 变形 。 
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图 4 法 向 过 载 时 间 历 程 


Fig.4 Normal overload time history 
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5 机 缀 B 时 刻 相对 A 时 刻 的 变形 
Fig.5 The deformation of the wing B-moment vs. A-moment 
表 1 机 辟 典 型 部 位 位 移 
Tab.1 Typical wing position displacement 


序号 民 向 位 置 弱 向 位 置 位 移 /m 
1 翼 根 前 缘 0. 033 
2 必 根 后 缘 0.026 
3 中 训 前 缘 0.169 
4 中 冯 后 缘 0.166 
5 辟 尖 前 缘 0.757 
6 辟 尖 后 缘 0.756 


根据 飞行 实测 的 数据 ,可 以 通过 两 种 方法 对 载 
荷 进行 静 气 弹 修正 :(D 是 根据 机 构 测 量 点 的 相对 变 
形 ,通过 搬 值 可 以 得 到 整个 机 忱 的 变形 ,进而 可 以 直 
接 得 到 因此 弹性 变形 而 引起 的 机 经 载 集 的 修正 量 ; 
G@ 是 通过 实测 的 飞行 参数 进行 载荷 计算 ,可 以 计算 
得 到 刚性 飞机 的 机 到 载荷 ,以 此 刚性 载荷 为 输入 , 按 
照 第 1 节 的 理论 可 以 计算 得 到 载体 修正 量 。 将 以 上 
两 种 途径 得 到 的 修正 量 进 行 对 比 , 可 以 根据 对 比 结 
果 对 方法 进行 验证 。 下 面 分 别 进行 叙述 。 


2.2 根据 实测 变形 修正 


用 B 时 刻 的 位 移 减 去 A 时 刻 的 位 移 ,得 到 在 1g 
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载 傈 基础 上 的 位 移 增 量 ,对 此 位 移 增 量 进行 插值 运 
算 , 可 以 得 到 机 可 每 个 气动 网 格 下 洗 控 制 点 处 的 位 
移 和 顺 气流 方向 的 斜率 。 插 值 方法 采用 薄板 样 条 插 
值 方 法 。 薄 板 样 条 持 值 函数 (有 文献 又 称 其 为 曲面 
样 条 函数 ) 是 径 向 基 函 数 的 一 种 , 径 向 基 函 数 是 一 类 
以 点 荆 到 节点 交 的 距离 d; = |x -x,| 为 自 变 量 的 
函数 。 利 用 径 向 基 函 数 , 将 试 函数 写 为 一 般 形 式 ” 


N 


u(xX) = > Pia pa (9) 
式 中 :9 = [o，w， py] 为 已 知 径 向 基 郴 数 ， 
本 研究 采用 薄板 样 条 函数 , 即 
bi(x) = dlogd? 
av] 为 待定 系数 。 根 据 定 


a = [ao wow 


解 条 件 
w(K) = W, (10) 
待定 系数 a, 可 用 已 知 的 节点 位 移 w = 
uy] 来 确定 。 在 本 方法 中 ,根据 测量 
点 的 坐标 和 各 个 个 测量 点 处 的 位 移 w 可 以 建立 一 个 
线性 代数 方程 组 ,求解 此 方程 即 可 得 到 待定 系数 
a; ,所 以 式 (9) 可 以 写 为 
u(x) =poa=po Hu=Nu (11) 
由 式 (11) ,根据 测量 点 的 位 移 可 以 插值 得 到 机 
翼 上 坐标 为 zx 的 任意 一 点 的 位 移 及 顺 气 流 方向 的 斜 
率 。 经 过 插值 计算 ,得 出 机 泌 实 测 位 移 和 气动 网 格 
控制 点 的 位 移 对 比 (图 6)。 然 后 计算 此 变形 引起 的 
气动 力 增 量 。 由 于 此 变形 的 存在 ,导致 飞机 在 同样 
的 飞行 参数 下 相对 刚性 飞机 总 升力 有 所 降低 ,为 满 
足 规范 规定 的 过 载 要求 ,需要 将 全 机 迎 角 统一 增 大 
某 一 角度 补偿 总 升力 的 降低 。 将 弹性 变形 引起 的 气 
动力 加 上 增 大 迎 角 后 增加 的 气动 力 即 可 得 到 最 后 的 
载荷 修正 量 。 
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图 6 机 可 实测 变形 和 气动 网 格 持 值 结 果 对 比 
Fig.6 Comparison between the measured deformation of 


wing and the result of aerodynamic mesh interpolation 


投稿 网 站 :http://cjam. xjtu. edu. en 微 信 公 众 号 :应 用 力学 学 报 


1268 应 用 力学 学 报 


2.3 根据 实测 习 行 参数 进行 修正 


在 试验 中 测量 得 到 了 飞机 J 迎 角 、 角 速度 和 舵 偏 
角 等 参数 ,由 这 些 参 数 计算 可 得 到 刚体 飞机 的 气动 
载荷 ,载荷 分 布 来 源 于 型 号 风 洞 试验 结果 ,机 翼 上 的 
气动 力 分 布 载荷 尺 是 速 压 dy 、 迎 角 a 、 侧 滑 角 B 、 舱 
偏 等 6 、 角 速度 p 等 飞行 参数 的 函数 , 即 
F = F(dg,a,B,6,,6,,p,q,r,…) (12) 
实际 型 号 中 ,都 有 足够 的 风 洞 试验 数据 ,根据 这 
些 数据 进行 必要 的 插值 ,就 可 以 得 到 的 机 可 上 某 特 
定 状 态 下 的 气动 力 载荷 分 布 。 同 时 根据 飞行 参数 
( 线 角 速度 、 角 速度 和 和 角 加 速度 ) 和 飞机 的 质量 分 
布 , 可 以 计算 得 到 飞机 的 惯性 载荷 ,将 机 豆 的 气动 载 
答 与 惯性 载 集 车 加 后 得 到 净 载 荷 ,以 此 载 丛 作为 输 
人 按照 第 1 节 的 方法 计算 即 可 得 到 根据 实测 飞行 参 
数 的 修正 结果 。 


2.4 两 种 途径 修正 结果 比较 

图 7 和 图 8 分别 给 出 了 机 小 前 力 和 弯 算 按照 实 
测 变形 的 修正 结果 和 实测 飞行 参数 修正 的 结果 。 图 
中 曲线 数值 分 别 代 表 剪 力 和 弯 矩 的 变化 量 , 横 坐 标 
表示 机 辟 展 向 相对 坐标 ,从 0 到 1 表示 从 必 根 到 
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展 向 相对 位 置 
图 8 弯 矩 各 站 位 修正 量 


Fig.8 Bending moment correction of every station 
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由 图 7 ~8 可 以 看 出 ,两 种 方法 的 修正 结果 趋势 
一 致 ,并 且 在 中 外 站 数值 上 也 比较 接近 。 而 在 靠近 
演 根 的 部 分 差别 稍微 增 大 ,这 是 由 于 在 实际 的 飞行 
中 ,飞机 很 难 做 到 真正 意义 上 的 稳定 俯仰 ,难免 带 有 
角 加 速度 和 结构 振动 ,在 笔 近 恤 根 部 分 因为 有 大 的 
集中 质量 (燃油 和 发 动机 ) ,导致 按照 稳定 俯仰 计算 
得 到 的 惯性 载 丛 与 实际 出 入 较 大 。 将 刚性 飞机 的 弯 
和 矩 与 两 种 修正 方法 得 到 的 弹性 飞机 的 这 矩 进行 对 
比 , 如 图 9 所 示 。 
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Fig.9 Comparison between the moment of rigid 
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aircraft and the modified moment 
由 图 9 可 以 看 出 ,两 种 修正 方法 得 到 的 弹性 飞 
机 的 弯 矩 符合 较 好 。 该 工 况 的 如 根 净 弯 和 矩 为 约 
2.3 x10”N. mm ,修正 量 大 约 为 3% ,修正 量 的 量 级 
也 较为 合理 。 


3 结 论 


本 研究 针对 大 型 飞机 机 四 气动 载 停 静 气 弹 修正 
问题 ,提出 了 一 种 计算 方法 ,并 通过 风 洞 试验 验证 了 
该 方法 的 准确 性 ,然后 将 该 方法 应 用 在 了 基于 试飞 
实测 数据 的 机 驾 载 荷 静 气 弹 修正 中 ,结论 如 下 。 

1) 针 对 不 同 刚度 的 机 总 风 洞 试 验 工 况 ,本 研究 
的 分 析 方 法 与 试验 结果 吻合 良好 ,证 明了 理论 分 析 
方法 是 准确 可 靠 的 。 该 方法 具有 物理 概念 清晰 、 计 
算 简单 和 精度 满足 工程 要 求 的 特点 ,适合 工程 应 用 。 

2) 完 成 了 基于 试飞 实测 数据 的 机 疲 载 集 静 气动 
弹性 修正 的 工作 ,以 两 方面 的 实测 数据 为 输入 ,通过 
两 种 途径 进行 了 修正 ,根据 机 可 变形 得 到 的 修正 结 
果 与 根据 实测 飞行 参数 得 到 的 修正 结果 符合 较 好 ， 
进一步 证 明了 理论 分 析 方 法 的 准确 性 。 本 方法 所 用 
到 的 机 辟 刚 度数 据 和 刚体 气动 力 数据 在 详细 设计 阶 
段 都 是 易于 获取 的 ,因此 本 方法 可 以 应 用 于 后 续 的 
型 号 研制 中 。 

3) 大 展 弦 比 后 护 机 中 的 气动 载 集 经 过 静 气 弹 修 
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正 以 后 ,可 以 明显 降低 全 根 的 弯 矩 ,本 研究 对 实测 试 
飞 中 的 特定 工 况 进行 计算 ,得 到 的 弯 和 矩 降低 量 是 
3% ,这 已 经 是 一 个 比较 可 观 的 数值 ,同时 该 工 况 并 
非 是 机 站 载荷 最 大 的 工 况 , 所 以 对 于 机 骂 的 设计 工 
况 载 荷 进 行 修正 后 , 弯 算 降低 的 量 级 会 更 大 ,这 对 于 
减轻 机 站 结构 重量 很 有 帮助 。 
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